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We derive the relativistic kinetic equation for Compton scattering of polarized radiation in strong 
magnetic field using the Bogolyubov method. The induced scattering and the Pauli exclusion 
principle are taken into account. The electron polarization is also considered in the general form of 
the kinetic equation. The special forms of the equation for the cases of the non-polarized electrons, 
the rarefied electron gas and the two polarization mode description of radiation are found. The 
derived equations are valid for any photon and electron energies and the magnetic field strength 
below about 10^^ G. These equations provide the basis for formulation of the equation for polarized 
radiation transport in atmospheres and magnetospheres of strongly magnetized neutron stars. 
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I. INTRODUCTION 
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1-^ ' Observations of the soft gamma-ray repeaters and anomalous X-ray pulsars showed that these objects can be 

^H, associated with the strongly magnetized neutron stars (NSs) with the magnetic field exceeding the Schwinger critical 
Q ' value of Bcr = m^c^/eh = 4.412 x 10^"^ G [1, 2]. This has revived the interest in theoretical studies of the interaction 
Vj , processes between radiation and matter in such fields [3]. 

C/3 ' Compton scattering is an important process shaping the radiation spectra of the NS atmospheres. Its properties 

, ^, ' in the magnetic field differs substantially from the case when the magnetic field is absent. Even the classical non- 

relativistic limit of the scattering cross section has a resonance at the energy related to the Lorentz frequency and is 

strongly dependent on the photon energy, polarization and the B-field strength [4, 5]. While the classical description 

^ . has been useful for understanding the approximate effects of energy, angle and polarization dependence of the cross 

'^ ' section in the magnetic field, it does not include the possibility of the electron excitation to a higher Landau state 

7%^ , corresponding to the resonances at higher harmonics, which required fully relativistic treatment. In the relativistic 

I~i ' regime the recoil of the electron is important and the natural line width of the cyclotron resonances depends on 

• [ the spin of the electron. The relativistic scattering cross section for the simplest case of ground-to-ground state 

^^ • scattering in the magnetic field was derived in [6] . These results were extended to a more general case of scattering to 

^^ [ arbitrary Landau states in [7, 8] and discussed further in [9] . The derived expressions have been applied to modeling 

the cyclotron line formation in accreting neutron star atmospheres, but only for the case of one-dimensional thermal 

electron distribution because of the complexity of the expressions [10-14]. When the incident photons propagate along 

the magnetic field, the resonance appears only at the fundamental frequency and scattering to the higher Landau 

levels can effectively be neglected. This allows to simplify the expressions for the relativistic cross sections and to 

; I , approximate them by analytical formulae [15]. 

Ci • The transport of photons through the atmosphere involves multiple scattering, which have to be considered either by 

the Monte Carlo methods or using the kinetic equations. The former approach was used for a qualitative study of the 

line formation process in Her X-1 [16, 17], but it becomes impractical for a large optical depth and when the induced 

scattering has to be accounted for, and therefore has a limited field of applications. In the cold plasma approximation, 

assuming the coherent scattering, the radiative transfer equation can be formulated as a set of coupled equations for 

two normal polarization modes [18]. The influence of the electron temperature on the radiation transport can be 

accounted by the Fokker-Planck approximation, for example, by modifying the Kompaneets equation [19] to allow 

for the resonances in the scattering cross section [20]. Such a treatment, however, does not account for the effects of 

the photon angular distribution and polarization. Photon polarization, however, influences the photon redistribution 
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over the energy [21, 22]. 

In a sufficiently strong magnetic field, owing to the large Faraday depolarization, the radiation can be described 
in terms of two polarization modes. Under certain conditions (depending on the field strength, photon energy and 
propagation direction), however, the vacuum resonance is accompanied by the phenomenon of mode collapse and 
the breakdown of Faraday depolarization [23-25] . In this case the two-mode description fails and instead the kinetic 
equations have to written in terms of the Stokes parameters or the coherency matrix. In the case when the induced 
scattering needs to be accounted for, the situation complicates further as there is no intuitive way to get such an 
equation. 

The aim of this paper is to derive from first principles a general kinetic equation for Compton scattering in any 
magnetic field accounting simultaneously for photon polarization in terms of the Stokes parameters, for the induced 
scattering and the Pauli exclusion principle for electrons. We use methods of quantum statistics and follow an 
approach similar to that used for derivation of the kinetic equation for Compton scattering without magnetic field 
[26]. The resulting equations are valid for any photon and electron energies, and for the magnetic field strength below 
about 10^® G. In the most general case, the electron polarization is also taken into account. We also consider several 
special cases and derive the kinetic equations when the electron gas is non-polarized and rarefied as well as when the 
radiation can be presented via two polarization modes. The derived equations provide the basis for construction of 
the models of radiation transport in atmospheres and magnetospheres of strongly magnetized neutron stars. 

II. DESCRIPTION OF THE ELECTRON AND PHOTON GASES 

We use the system of units where h — c ~ rric ~ 1. We assume that the magnetic field is locally homogeneous, 
which is justified, because the scales of changes of the B-field are orders of magnitude larger than the microscopic 
magnetic scale for conditions in atmospheres of NS and even the geometrical depth of the atmosphere. The magnetic 
field is described by the dimensionless parameter b — B/Bcr- We choose the reference frame in any space-point so 
that the z-axis coincides with the magnetic field direction. The following assumptions about the time scales are used: 

1. The typical time scale on which the distribution function changes (for electrons and photons) is much larger 
than the typical time scale between the interactions. 

2. The plasma is sufficiently rarified, so that we can use a generalization of the Bogolyubov method for the case of 
quantum statistics to derive the kinetic equation. 

3. The typical time scale of a single interaction is much smaller than the typical time scale between the interactions. 

A. Descriptions of single particles 

The electron states are described by the wave- functions ^no-dl, Y, Z). Its arguments are the space-time coordinates, 
the momentum projection Z on the direction of the magnetic field, Y describing location of the center of electron 
gyro-orbit (its y-coordinate), the Landau level n, and the spin projection a on the magnetic field direction (cr = ±1 
in h/2 units). Dimensionless energy of an electron in this case, 

Rn{Z) = ^/l + Z^ + 2bn, (1) 

is independent of Y. The energy levels are degenerate with the spin projection on the magnetic field direction, except 
for the ground Landau level with n — 0, where a can have only one value —1. The full electron wave function is 
presented through the partial solutions of the Dirac equation for the electron in the magnetic field: 

where 6„o- are coefficients. 

The photon state is described by four parameters: the wavenumber k, the two angles 9 and ip, which define the 
direction of the photon momentum, and the polarization state s = 1,2. The 3-dimensional photon momentum can be 
represented as 

k = {kx, ky,kz) — fc(sin0 cos(p, sin sin (p, cos 6). (3) 



The corresponding photon 4-momentum is fc = {fc, fc}, fc = |fc|. Photon polarization is described by the polarization 
basis. It consists of two unit vectors, which are orthogonal to the photon momentuni k: 

ei ~ {sinip, — cos(p,0), 62 = (cos^cosi^,cos^sin(^, — sin^). (4) 

The 4- vector potential can be defined as: 

Asiz) = e,e-^--, e,:={0,ej, s^l,2. (5) 

We note that the photons are described in the same manner as in the case when the magnetic field is absent. We 
assume that the dispersion relation for the photons in magnetized vacuum does not differ from the dispersion relation 
in the case when the magnetic field is absent. This approximation constrains the strength of the field and energies of 
photons. For estimations one needs to know vacuum dielectric tensor and the inverse permeability tensor for the case 
of magnetized vacuum [27, 28]. It is known that the indices of refraction differ from unity by more than 10% only for 
the fields with strength b > 300 [29]. This restricts application of the developed formalism to B ^ 10^^ G. 

B. Description of the particle ensembles 

1. Wave functions 

We describe particle ensembles by density matrix using the rules of quantum statistics. Let us define the wave 
functions for the case of limited number of particles. These functions will be used for construction of the density 
matrix. The wave function for a limited number of particles with defined characteristics of each of them, can be found 
from the vacuum wave function by applying the operators of creation and annihilation. Let arg) (k) and au\ (k) be the 
creation and annihilation operators of a photon in the state with polarization s and 3-momentum k. According to 
the methods of second quantization [30] , these operators satisfy the relation 

a(,) (fc)a(,,) (fc') - a(,,) (fc')a(.) (k) = kS{k - k')S'J . (6) 

Let 5jjg.(Y, Z) and bna(Y, Z) be creation and annihilation operators of an electron on the Landau level n in polar- 
ization state (7 with momentum projections Y and Z. These operators satisfy the following relation 

bnAY,Z)bl^,{Y\Z')+bl^,iY',Z')bUY,Z) = R,,{Z)SiY-Y')5{Z-Z')S:i'5^J. (7) 

The system of N photons with fixed parameters {si, fci; ...; sn, k^} is described by the wave function 

*si...s„(fci, •••, few) = ^=a(«i)(fci)...a(^„)(fcAr)*o, (8) 

where ^Pq is the vacuum wave function of the photon gas. Analogously, the system of N electrons with fixed parameters 
{rii, CTi, Yi, Zi; ...; nj^, un, Yn, Z^} is described by the wave function 

K\::.ZiYuZu...,YN,ZN) - -^bl,^iY,,z,)...bi^^jYN,ZN)^„. (9) 

where $0 is the vacuum wave function of the electron-positron gas. The wave function for arbitrary state of particles 
can be presented as a sum of the wave functions with fixed particle parameters. For example, the state of N photons 
is described by the function 

f ^ dk 

*^= / '[[j^c,,,„s^{ki,...,kN)-^s^...s^{ki,...,kN), (10) 

where Csi...sjv(^i: •••: ^^n) are the weight coefficients. The wave function for an arbitrary state of N electrons can be 
represented as 

The wave function for state with A^ photons and A^+ electrons can be written as 

T / 1 — r ^"^i 1 — r G.Y jQ.^ j o"! , , .(Tjv 1 , . ., 

^N,N+ ^ / 11 ~T~ 11 "^ 7yTCni...nN+.si...SN(Xl^ ^1^ • • • , ^A^+ , '2^A^+ , ^1 , •••,%) 

i—1 j — 1 3^3' 

X ^,^,„s^{ki,...,kN)'K\:'ll\{Yi,Zi,...,YN^,ZN^). (12) 



2. Density matrix 

The density matrix is defined as an averaged dyad product of the state vector with its conjugate. For the system 
consisting of A^ photons and iV_^ electrons it can be written in the form 

1 /"T^dfc.dfe' r:^ dY,dZ, dY,'dZ' 



p.,..=^7r^/ni^9^ n"""^^ 






_ (T-, ...a-f. 



The expressions in the triangle brackets are the elements of the density matrix kernel. 

3. Algebra of the density matrix kernels 

From now on we will operate only with density matrix kernels. All the equations and the final results are written 
through these kernels. It is easy to make a transformation from the simplest kernels to the distribution functions or 
to the coherency matrix. Let us write the density matrix kernel for the system of N photons and N^ electrons 



Psi...sjv,o'l---CTjv^,ni...njv_|_ I , , 



Yi...Yn^ Zi...Zm+ 



C'„/ „' „' „/ (^1, -Z^l, •••, yAT, -Z^AT, , fcl, •••, fcAr)Cni...nN, ,si...sjv(^li ^1: •••: ^JV+^AT+i fcl, •••, ^^Af) ) • (14) 



n' ...n' ,s'...s' V"- i ' "i' •••' " iV+"iV+ J '"iJ •••! '"iV y'-"i---"N_|_,»l---sjv v-" J-i "J-i ■■■^ -'iv+^jv+j 

Further let us write some useful relations for the kernels. For the sake of simplicity we consider only the photon gas. 
These relations can be generalized trivially to the case of the electron-photon gas. The kernel for the system of TV 
photons can be written through the density matrix: 

Hereinafter we call it the A^-particle kernel. It is normalized to unity: 

nt<:::s(:::t)-p(^)^i- («) 

For any m < N the m-particle kernel can be calculated as 



P'^-'- W-hn (N-m)\ . ii h ^^---^-+--" Ui-femfcm+i-fcA^y ■ ^ ^ 

The 1-particle kernel can be expressed through the A^-particle kernel as 

k) {N-iy.J l\ h ^' ^^■■"" I k k,...kj, ) ^'"^ 

It is normalized to the total number of the particles: 

d^fc „ fk 



The diagonal elements of 1-particle kernel compose the coherency matrix in the case of photon gas. 



4- Transformation from the simplest kernel to the distribution functions 

The transformation from 1-particle density matrix to the distribution function in the case of electrons in the case 
of field-free space can be made using the Wigner function. The Wigner function is defined as 

p{p,r)= dvcxp (ip ■ v) ps {r+ v/2,r— v/2) , (20) 

where pg (r+ v/2,r— v/2) is a 1-particle density matrix in the coordinate representation. The momentum and 
coordinate representations are connected through the Fourier transforms: 

, 1 f dp dp' I /\\ fp'\ 

Ps{r,r) = -——r^^==cxp{-t{p-r-p-r))pl . (21) 

Then one can rewrite the Wigner function using the density matrix in the momentum representation: 



^(^' ^^ = (^ / ^S^ '^"p ^-'^""^ - ^^^ • ^) ' (^ - 



^^+^^VrM. (22) 



2 J \Pi 
The inverse transformation from the Wigner function to the density matrix in momentum representation reads 



p(^;)=^|^/dpdrexp(»(p,-K)-r')-5(p-^l±^)p(p,r). (23) 

The time scale of the electron-photon interaction is much smaller than the time scale of noticeable changes of the 
distribution functions. Therefore, in the last equation one can assume that the Wigner function does not depend on 
the space variables. In this case, the integration could be made easily and we can write 

p(pj)=m'5(W-Pi)p(Pi). (24) 

Then one can convert the 1-particle density matrix kernel to the distribution function in the case of spinless particles 
or to the coherency matrices in the case of particles with non-zero spin. In the last case a trivial generalization is 
used: 

Pli (^1) ^PoiS{p[-p,)pll{p,). (25) 

Equations (24) and (25) can be written immediately from the physical meaning of the density matrix and the assump- 
tion that the time scale of the interaction and the time between interactions are mush smaller than the time scale of 
noticeable changes of the distribution function. 

One can also write similar relation for the case of electrons in the external magnetic field. If we consider a non- 
interacting electron in the B-field not accounting for cyclotron radiation (which should be described by another kinetic 
equation), then the electron should conserve its z-momcntum and the Landau level. This means that the kernel should 
be diagonal over both Z and n, because it is not possible to have mixed states corresponding to different values of the 
2;-projection of momentum or the Landau level. Non-diagonal elements in the kernel can appear only if one accounts 
for interactions between particles, but because of the smallness of the interaction time scale the kernel should be 
diagonal over z-projection of momentum and the Landau levels. In this case the relation will have the following form 
(a detailed derivation is given in Appendix A): 

pr (^^1') =i?„(^)<5::'<5(y'-r)(5(z'-z)p-l(y,z), (26) 

where a and a' describe the electron spin-states, n and n' arc the Landau levels, Z and Z' arc the momentum 
projections and Rn{Z) is the electron energy given by equation (1). 

A transformation from the 1-particle density matrix to the distribution function in momentum space is trivial, but 
one must again assume that the typical time scale of changes of the distribution function is much larger than the 
typical time scales of interaction between the particles. 



C. Description of the interaction 

1. Description of the single interaction 

Let US mark parameters of the particles before the interaction with the subscript " i" and particles after interaction 
with the subscript " f ' . There are three conservation laws for Compton scattering in the magnetic field. They are 
the energy conservation, the conservation of the momentum along the magnetic field and the conservation of the 
transversal momentum: 

-Ri + fci = -Rf + fcf , Zi + ki cos 6i = Zi + fcf cos 0f , Fj + fcj sin 6i sin (fi — Yf + k{ sin Of sin ipf. (27) 

Let us use special designation for product of (5-functions which are describing these conservation laws: 

6{n,Y,Z,k\ n',Y',Z',k') = 
= S{RniZ) + k- Rn' {Z') - k')5{Z + kcos9- Z' -k' cos 9')6{Y + fc sin sin (p - Y' - k' sin 6' sin (p') . (28) 

A single interaction can be described by the S-matrix. The elements of the S'-matrix can be calculated using methods 
of quantum electrodynamics. In the simplest case the elements of the S'-matrix can be obtained using second-order 
perturbation theory. In this case Compton scattering can be represented by two Feynman diagrams with two vertices 
in both of them and one can write an expression for the S'-matrix elements: 

^fi = ^Ama j &\^&\2^f{L2) { [i4{L2)\ G{l2,Li) h^Li)] + [iMl2)\ G{l2.Li) [zMiLi)] } *i(z:i), (29) 

where 7 A is the Dirac inner product of a 4- vector and 7-matrix, and G{z2tLi) is a relativistic electronic propagator in 
the presence of a constant magnetic field, ^i(r) and ^f(r) are the electron wave- functions in coordinate representation, 
and a = e^ is the fine-structure constant. The S'-matrix elements and the cross-sections for Compton scattering in 
magnetic field contain resonances which have to be regularized [31]. The calculations are not trivial and have been 
performed only in special cases [6-8, 15]. 

2. Evolution of the density matrix 

The evolution of the density matrix can be described by equation: 

^^ - H{t)p{t) p{t)H{t), (30) 

where the Hamiltonian is 

H{t) = -e / dr^{r)2_A{r)ij{r). (31) 

Equation (30) is written here in non-covariant form, but it can be transformed to the explicitly covariant form using 
the Tomonaga-Schwinger equation [30]. It means that the form of the equation is covariant for the longitudinal 
Lorentz transformations (along the magnetic field direction). The solution of equation (30) can be presented by the 
operator of evolution C/(x, y): 

t 
J ^^'^^' " '"^^^ " ''^*°^ " ^^*' *°^'°^^°^ " Piio)U{t, to). (32) 

to 

On the other hand the operator U{t,to) can be represented in the following form: 

t t 

U{tM) = /di' /dt" j dr' I Ar"S{r!,r!') = /dV f d*r"S{r' ,r:'), (33) 

to to V V V V 

where V = [io , ^] x ^ is the volume in Minkowski space and 

dYdZdY'dZ' dkdk'. _„, ,_ ,.„ ,,..r^>„'.' f k'Y' Z' 



^, , ,,, e^ f dYdZdY'dZ' dkdU . , ,,, , , ,,,, ... . r„'^'.' / 'fe' 
5(r',r") = ^J-^j -^ ____fo1^,^,(z')6„.(Z)a(,,)(fc')a(.)(fc)AA;j ^ ^ ^ 



Y Z 



(34) 



The space integral of N can be represented through the elements of the scattering M-matrix: 



dV / dV'AA: 

V V 



4^// \f7l<JS 

n a s 



k' Y' Z' 
k Y Z 



{2n fS{n,Y,Z,k\n',Y\Z',k')M^/ 



n'Y'Z' 
nYZ 



(35) 



then U{t,to) can be rewritten in the following form 

dYdZdY'dZ' dkdk' . 



ZTT 



R 



R' 



k k 



rblAZ')hnAZ)-a(s'){i^)a(s){k)M^J/ 



n'Y'Z' 
nYZ 



(36) 



Let us assume that the typical time scale of the density matrix changes is much larger than the typical time scales of 
a single interaction. In that case changes in the distribution on a macroscopically small times scale can be represented 
througth the S'-matrix because M-matrix can be considered as the scattering S'-matrix divided by the fine-structure 
constant: 






= Mfi- 



Sj\ 



(37) 



and the time interval [ioi^o + i] is considered as a macroscopically small time. Equations (32) and (36) determine the 
solution formulated through the elements of the scattering matrix. We reformulate these equations below in terms of 
the kernels of the density matrix. 



III. DERIVATION OF THE KINETIC EQUATION FOR THE PHOTON GAS 

A. Methodology of the kinetic equation derivation 

1. Summary of our assumptions and the Bogolyubov method 

We derive kinetic equation using a generalization of the Bogolyubov method (for the case of quantum statistics). At 
the first step we formulate Liouville's theorem in terms of the kernels of density matrix. One can derive the equations 
for kernels of different orders (1-particle, 2-particle and other) by integrating over the parameters of different numbers 
of particles. We use this method to obtain the system of kinetic equations. If the full ensemble contains N particles, 
the system of equations contains N equations. In the case of the rarefied gas one can use only a few first equations 
from the Bogolyubov hierarchy. The criterion of rarefaction can be formulated through the "gaseous parameter" , 
which depends on the concentration of the particles and the cross sections of their interaction: 

agas = V^inenp^)^/^ < 1, (38) 

where ctt is the Thomson cross section, n^ and nph are the electron and the photon concentrations, correspondingly. 
According to the principle of weakening of correlations, which is satisfied for sufficiently rarefied gases, the correlations 
are accounted for only in the equation for the 1-particle matrix via kernel by entering the right-hand side (rhs) of the 
aforementioned equation. This kernel is assumed to characterize the electron and photon states after the interaction. 
It can be represented via the same kernel before the interaction and the correlation function. To derive the kinetic 
equation for the typical conditions in the neutron star atmospheres, it is enough to use only the first and the second 
equation from the Bogolyubov hierarchy. 



2. Formulation of Ltouville 's theorem and the equations of Bogolyubov hierarchy 

We use the following notations: R = Rn{Z), R' = Rn>{Z'), etc. There are N photons and iV+ electrons in the 
system. The equation, describing the change of {N + A^+)-particle kernel during macroscopically small time Tq is 
written as 

2 
dYdZdY'dZ' dkdk' 
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(39) 



This equation is a formulation of Liouvillc's theorem in terms of the density matrix kernels. It can be used to find 
the equation describing the evolution of 1-particle photon kernels on the time interval [—To/2, To/2], which is the first 
equation of the Bogolyubov hierarchy. 

We integrate equation (39) over the parameters of {N — 1) photons and iV+ electrons: 
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where the kernels under the integral in the rhs of the equation correspond to time t = —To/2 and where we used 
special designations for kernels describing the system of A^ photons and iV+ electrons: 



psa ■ ^' 

'' ^ kd 



i,J, 
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where J ~ (s, k) and G ~ [a, n, Y, Z) are parameters of photons and electrons respectively. Let us denote the terms 

N 

under the sum ^ by S. One can transform it: 

i=2 



1 



^^ f^K^'^r^^^Sin,Y,Z,k\n',r,Z',k') 



{N- iy.N+\ 2tt J R R' k k 

'n'Y' Z' 



n Y Z 



k'\ ^ f4L/dk,dk'^ ^+ 



i+ = \'' i=2 ^ * ^ ' i+ = \ \ *+ 



X hi.„. (r', Z')6„,(y, Z)-a,, ik')asik)^C 

( , " ' / k' 

X }^5:,^S^^,^5iY' -YpSiZ' - Z',^)J2SI^S{k' - k,)P!l (^ ^| 

- SlJ^JiY Y,^)SiZ Z,^)J2stAk-h)Ri i ^1 



E/ni^^ini'''^'^^^''^''^^^ 



K 



i,i+,J, G 



i=2 



N+ ,. N 



i+ = l 4=2 



{N - l)\N+\ 27r 

N+ 



i,i+,J', G" 

/dkdk' 
dYdZdY'dZ'- — —S (n, Y,Z,k\ n' , Y\ Z' , k') 
k k' ' 



y /T-r/'d^d^\ 1^ {dY,^dZ,^dYl^dZ[^ 



N 



R. 



■1+ 



i?' 



.4 = 2 



ki 



i, t+, , 



j\g)-j:phI: 



i=2 



ki 



1,1+,. ]\G' 



Creation and annihilation operators with the elements of scattering matrix were placed under the integral in the last 
transformation. Indices in the round brackets indicate the positions of pairs of parameters {s-s'), {n-n'), {Y-Y')^ 
{Z-Z'), {(J-a') and (fc-fc'). One notices that all the terms in the sum S cancel out, giving S = 0. The remaining part 
of the equation can be rewritten as 



Psi 



; (K 



'\h 



,^<,.,. , n'Y' Z' 
'"" ' n Y Z 



u 



^\-Ps 



Ta\ a f dYdZ dY'dZ' dkdk' , , , , , ,n 

^ i— I ^ :r^ -^j-5{n,Y,Z,k\n',Y',Z',k') 



271 R R' k k' 



k 



xs' xnJ iJ \ „san 



s:,s(k'-k[}pr^;.,.i^^ 



Z) -'■■■''"- "'^"i'-'-Ci' 



YZ 
Y'Z' 



(41) 



Thus, we have obtained the first equation of the Bogolyubov hierarchy describing the evolution of 1-particle density 
matrix kernel through the 2-particle density matrix kernel. 

Let us now obtain the second equation of the Bogolyubov hierarchy for the 2-particle kernels. We proceed with the 
integration and the summation over the parameters of (N — 1) photons and (N^ — 1) electrons in equation (39): 






To\ _ s[cr{n[ (kiYlZ[ 

2 J ^^^"^"^ Ui^i^i 



To 



1 



{N - l)\{N+-l)\ 27r 



dYdZ dY'dZ' dkdk^^ ^^^ ^^ ^ ^ I ^,^ ^,^ ^,^ ^^ ^,^^^^^^,^ ^,^^^^^^^ Z)-aA1^)a.{k)M:'/ ('^' ^' ^' 



R R' k k' 



^ '- '--'^ ""' ' dY,^dZ,^dYl^dZ[^ 



ni^fin 



R 



m 






k 



St, 



[s{k' k[) ^ </;';/(K,, - Y')s{z,^ z') p-:^ [^ ^^ 



s:i,s:,st,5{r-Yi)s{z'-z[)s{k'-k[) p:i^^ 



l,i+,J,G 



1,1, J, G 



i+=2 



N 



5ifJ,5{Y' Yl)5{Z' Z[)J2si5{k-h) P:^ 



i=2 



i, 1, J, G 



N N+ 



i=2 i+=2 



1,1+ 



i+, J, G 



SI'5';,'S^'6{k-ki)5{Y-Yi)5{Z-Zi) R^j 



u 



1,1,J', G" 



St^d{k-h) J2 S"S:^d{Y-Y,AS(Z-Z,A Rl) : 



i+=2 



1,*+,J', G' 



N 



■ 5:,-5%-5{Y -Y,)5{Z - Z,)Y,KAk- h) Rl\ 



i,l,J',G 



N N+ 



E E '^sV«./^./(^ - Y^+)SiZ - Z.^Sik-h) Rll 



i=2 i+=2 



i,i+,J', G' 



(42) 



AT N+ 



The terms with a double sum 2_. /, cancel each other. Other terms are transformed into the form containing 2- 

i=2 i+=2 

and 3-particlc kernels. As a result, we obtain the equation for the 2-particle density matrix kernel: 





s\ g\ n\ I "'1 


Y{Z\ 


To 


\ s^a\n\ I «! 


y;z\ 




T 


A 




Ps,.,n, y ^^ 


YiZ^ 


2 


) ' ^^^'^^"^ ^v fci 


YiZ^ 


-^J 




a f dYdZ dY'dZ' dfc dfc' , , 

— / In, 

2tt J R R' k k' ^ ' 


, , f -n' Y' 7J 

Y,Z,k\n',Y\Z',U)M^„; ( ' X 

\n 1 A 


k' 
k 


X 


8i,^si.f;^6{Y' -Y[)6{z' 


-Z[)S{k'-k[)pl-\^f^ 
\ "'1 


YZ \ 
YiZj 




- 


^K.KA.A^ 


- YxA 


Z- 


-Z,)Sik-h)pl 


a'n' 

\ 


k' 


1 


Y{Z{\ 
Y'Z' ) 
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YZ 

YiZi 



kik' 

Y{Z[YZ 
YiZiY'Z' 



-51X^5{Y-Y,)5(Z-Z,) pI\ 



St S(7 n 

s' at n'-. 



kiU 



SlAk-k,)pXi:i 



St (Tt rii an 
a' n' 



YlZ[YZ 
YiZiY'Z' 



Y'Z' 
YIZ[ 



(43) 



Thus, we have derived the equations for kernels of 1-particle density matrix of photons (41) and for the 2-particle 
density matrix (43), which includes both the photon and the electron parameters. 

B. Completion of the derivation 

1. Expression for the 2-particle kernels through the 1-particle kernels 

Equations (41) and (43) are the first two equations of the Bogolyubov hierarchy. The hierarchy can be continued, but 
using the principle of weakening of correlations, we have stopped at the first two equations. To obtain an equation 
describing the evolution of the 1-particle kernel through the 1-particle kernels, one must use these two equations. 
We use the "molecular chaos" approximation, according to which there is no correlation between the distribution 
functions of photons and electrons before interaction. This approximation works better in the case, when the typical 
time between the interactions is much larger than the typical time of an interaction. The independence of the photons 
and electrons distributions can be expressed through the following equation: 



Pa 



Y'Z' 
YZ 



= Pa 



Y'Z' 
YZ 



(44) 



The 2-particle kernels of the photon gas and the 2-particle kernels of the electron gas are presented through 1-particle 
kernels according to the properties of symmetry and anti-symmetry of bosonic and fermionic wave functions: 



PsiS2 



cr^njCTjrij / Y1Z1Y2Z2 

P<TiniCr2n2 [ ^r ry ^r ry 

Y\ZjxY2^2 



k\k2 



= Ps\ 



- rt"!"! 

— Pairii 



Y{Z', 
YiZi 



-Pst 



„'^2"2 
PO"2n2 



YoZ-2 



„0'2"2 

PcJini 



k[ 
h2 

^2^2 

YiZi 



P<T2"2 



^2^2 



(45) 
(46) 



These equations become more accurate, when photons and electrons gases are sufficiently rarefied. Transformations of 
the 3-particle kernels in equations (41) and (43) are simple because there are no pure photon or pure electron kernels 
among them and one can rewrite them easily through the 1- and 2-particle kernels. Then we use equations (44) and 
(45) to complete the transformation. 



2. Closure of the Bogolyubov hierarchy 



Substituting equation (43) to equation (41), we get: 



Ps{ 



fci 






To 
2 






dYAZdY'dZ' dkdU 



R 



R' 



k k' 



S{n,Y,Z,k\n',Y',Z',k')M^'/ 



n'Y'Z' 
nYZ 



k'5i,{u-u,)pii:,^, 



fcl 



^fj-fc'^^^cfc-fciv:?;;- 



fc' 



YZ 

Y'Z' 



E 



AYdZ AY'AZ' dY"dZ" dY"'dZ"' dkdk' dk" dk" 



(27r)2 ^^ J R R' R" R'" k k' k" k" 

n,n' .n" .n'" 

X 5 (n, Y,Z,k\ n', Y' , Z', U) S (n", Y" , Z" , k" \ n'" , Y'" , Z'" , k" 



X m;/ 



n'Y'Z' 
nYZ 



M ,, f, 

-"^-*cr"s" 



n"'Y"'Z"' 
n"Y"Z" 



k'" 
k" 



xlsiik'-k',) 



'IJI'TJI'IXS n (J 



k'k"'R"'6 



{k-k"')6{pyz-p'{^z)pk:'r^ 



ki 



Y"Z" 
Y'Z' 



11 



k'k"R"5i[::f{k, k")s{p'yz p'{-z)pt^'."'n"' 



k'R"'5-„ - (p-^ - PYz)pZs!^'a'n' 






fcfc" 

fclfe'" 



Y'Z' 



YZ 

\'" 7,'" 



k'k"'5t {k-k'")pl,:r^% 



YZV'Z" 
Y'Z'Y"'Z"' 



YZY"Z'' 



Y'Z'Y"'Z"') ^'^"^"n'^' ^'P'^Z PYz)pTis""an [^^g 



kk" 



Y"'Z" 
YZ 



-kk"'B"'5i":i"^"'{U"-u,)5{py^-p'i^,)pi:::y ( ^, 



5l^{k-k,) 

+ kk"R"5i!{U - U')5^n'ar"''"{p'yz PvzV^Z^ 

1 1 1 II 



Y"Z" 
Y'Z' 



k[ 
k'" 



YZ 
Y"'Z'' 



kR"'5iy"'{py,-p'i^M^-,:' [ li^iii 



Y"Z' 
Y'Z' 



kR"6-,J, {p'y, P^z)pX'' 



s' s'" an 






YZ 



-kk"'s:::'{k[-k"')p: 



/'/\ ^s' a n 



k' 



YZY"Z" 
Y'Z'Y"'Z"' 



+ kk"^" (k' - k"]n^ 



k[ 



YZY"Z" 
Y'Z'Y"'Z"' 



, (47) 



where we introduced a symbol for the product of several Kronecker's deltas: (5^^' "q^ = 11^=1 '^a;' ^ symbol for 
the product of Kronecker's deltas and a (5-function: 5a\"'al,i'^) = '^Qi"'Q"'^('^)i ^^.^ Pyz — (Oi ^i ^) is the electron 
momentum. 

Using equations (44), (45) and (46), one can transform 2- and 3-particle kernels in the rhs of the equation to the 
1-particle kernels. Then using equation (25), the equation can be represented in terms of the coherency matrix. After 
some algebra we get: 



Psl 



; fk[ 



ki 



To 






2 



9^ 2^ 



27r 



a 
4^ 



E 



drdz dr'dz' dfc dfe' , ., , , , .n . 

drdz dy'dz' dy"dz" dY"'dZ"' dkdk' dk" dfc'" 



n.7i' ,n" .71 



R R' R" 

X S (n, Y,Z,k\ n', Y' , Z', U) M^' 



R'" k k' k" k" 
n'Y'Z' 



nYZ 



X 5 (n", Y", Z", k" I n'", F'", Z"\ k'") M^',','/," 



n"'Y"'Z"' 



n"Y"Z" 
X {Bi +82+83+8^ + 65+ 8e} , 



k'" 
k" 



(48) 



where 



Ai = m;; 



n'Y'Z' 
nYZ 



k' 



R'Sl {pyz - p'yz)p:'n'iPYz) k' 5{U - k[)5{k-k,) (k,Sl} pl^{k,) - kSt^pl}{k')) ,(49) 



8, = k'k"'6{k" - k"')6f{pyz p'^zK'iPYz - p'^z)P's'"ik"') (R'R"'6fp^;MYz) RR"5^j: pC{Pyz) 



X (k5{U, - k')St^{k-h)pl}{k') - k,S{h - k)Sl[{k' - k[)pl^{h)) , 

82 - R'S:'{pyz p'y,)p%,Ap'Yz)R"'6i'"{PYZ PYz)P:''"n"'iPYz) 

X r5^< (fc' - k[) {k'k"k"'5{[ (fei - h!')5{k - U")pI,„ {k'") - kik'k"'Sf' {k - k"')S{k" - ki)pl" (fei 
+Sl^{k-ki) (kk'k"S'J,"{k[ - k"')5{k' - k")pl','{k') - kk"k"'dl"{k' - k")5{k[ - k"')p%{k"')] 
8s - St^ik- h)kk'R'R"'5{k' - k")S{k[ - k"')S:::{p'yz - p'^z)Sf'{PYZ - PYz)pi: {^)p%Upyz) 
X 'k"'5fp%{U") + y"5i;' (6f p:,„^,4Pyz) 



(50) 



(51) 
(52) 



8, - -k,k'k"'R'R"'5:,{k' - k[)S{k^ k"')6{k, - U')5l {pyz - p'^z)K' (p'yz - Pyz)pI, {ki)p%in'{p'Yz) 



B, = k'k"'5{k U")5{k, U')5i. {k' k[)Sf (pyz p'{;z)S:; (p'yz p'^z)pl,„ {k'") 



X RR"k^5:, pl^ (fci)p-„ (pyz) + k"R"'5l^ pl,„^,„ (p^'^) [B!'8l, - R'p:,APyz) 
B, = -kk"'6{k[ k"')S{k' k")dt^{k-h)df{pyz p'^^zK'iPYZ Pyz)p%{I^ 

X 'RR"k'5%:pi:{U)pC{PYz) + k"R"'5i:p%,„^,„{p'^z) (^"5%: - r'p%:apyz) 
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(53) 

(54) 
(55) 



3. Simplification of equation (48) 



The presence of the J-functions under the integrals in equation (48) allows us to reduce a number of integrations. 
Let us define two singular measures pi and p2: 



drdZdr'dZ'dfcdfc' , ,, , , , ,„ 

dpi = — — —5{n,Y,Z,k\ n\Y',Z\U) , 



R 



R' k k' 



^^ R R' R" R'" k k' k" k'" " 

S (n, Y,Z,k\ n', Y' , Z' , k') 6 (n" , Y'\ Z", fc" | n'", Y'" , Z'" , k'") 



The terms from the rhs of equation (48) can be written as: 
V- /■ .NT- f dYdZ ,, (nYZ 



Ejd^^^i = E/ 



R "'"'' \nYZ 
dYdZ dY'dZ' dfc 



p-,„(z) 5t,^pt^{ki)-st,p:}{h 



R R' 



-S(R' - R)M^,, 



nYZ 
n'Y'Z' 



k\ ,, n'Y'Z' 
k ^' \ nYZ 



(56) 



(57) 



(58) 



p\,„{k) 



[sfpUAz') - €pC{z)] [siykK) St[pt^{h) 



(59) 



E j^P.B, = Y.j 

x(5ilpl,„{k)-5fpi[{ki 



dYdZ dY'dZ' dk 



R R' 



5(ki-k)p:,Az)p:,„„,(z') 



m:, 



1 ^ 


(nYZ 


M'^ " 






\nYZ 



k\ ,,a"'s"' f n'Y'Z' 
kj " ' \ n'Y'Z' 



,; fnYZ 
\nYZ 



fci\ '"'" f n'Y'Z' 
k / ""'" \ n'Y'Z' 



-S^'p^4ki) 



(60) 



X M^ 



E 

n,...,n 

n'Y'Z' 



j dp^B^^Y^j 



dYdZ dY'dZ' dk 



R R' 



6{n,Y,Z,ki I n',Y',Z',k) 



nYZ 






7i,n 

nYZ 
n'Y'Z' 



''^^ pi:{k)p:i,(z') [sfp%{h) + str [sf - p^^,„jz))] , (6i) 



n,...,n 



Ejd^^^^ = -E/ 



dYdZ dY'dZ' dk 



R R' 



S{n,Y,Z,k\ n',Y',Z',ki) 



""^^"^ ' nYZ 



k / " ' \ n'Y'Z' 



ki 



pl^{ki)p%.AZ') S: p%4k)+St [S: -p%„.AZ)\ , (62) 



Ejd^.^a = E/ 



drdz dy'dz' dfc 

~R R/~ 



S{n,Y,Z,k\ n',Y',Z',ki) 



xM-^'-'^'^' 
' nYZ 



ki 






nYZ 
n'Y'Z' 



f^]pi'"{k) [s:':pi';{h)pC{z)+si';p:,„jz) (^s^: - p:i,{z'))] m 
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^ Jd^2fi6 = -^| 



&YdZ AY'AZ' Ak 



xM^ 



nYZ 
n'Y'Z' 



n,...,n 

k 



h 






n'Y'Z' 
nYZ 



fci 



R 



P%{ki 



R' 



5{n,Y,Z,k\ n',Y',Z',ki 



C pl, {k)p%AZ') + St, P%.,.,AZ') [5%, pZ,^{Z) 



(64) 



We use the rhs of equations (58)-(64) in the rhs of the final kinetic equation. One can notice that the rhs of expression 
(59) vanishes after summation over the electrons spin states and the photons polarizations. 



4- Transformation of the left-hand size of equation (48) 

It is necessary to rewrite the Ihs of equation (48) in terms of the distribution function. Then one can rewrite Ihs 
through the differential operator. The later transformation is similar to the one, which is made in the quantum field 
theory for transition from the limited space-box to the infinite space. Finally we get: 



Psi 



' / k' 

1 I ^1 



Psi 



' /k' 



To\ d s[ (K 

-y)=^°d^^-Ui 



: rofci<5(fc'i - k^)-pl\{k,) = k,6{U,~h) — —pl\{k,), (65) 
At c At 



where in the later transition the relation 5{ki — k'l) — cTq/{2tt) is used. Then one can restore the dependence of 
the photon matrix on time and space coordinates. After rewriting the derivative over the line of sight as the full 
derivative, the Ihs of the equation takes the covariant form 



2tt Ar ' / ki d \ ' ' 

ki6{k:,~k,)——ptl{h) ^^ 2^,5(fci - fc'i) ( ^- + fcij ptlik,,n,t) = 2TTSih- k:,)k,Vptiiki,n,t). 



(66) 



IV. DIFFERENT FORMS OF THE KINETIC EQUATION 
A. Kinetic equation for coherency matrix 

1. The general form of the kinetic equation 

Now one can write the final form of the kinetic equation. In the most general case, we formulate it for the coherency 
matrix, where the polarization of electrons is taken into account (i.e. for the situation when there can be non-trivial 
spin-distribution of the electron gas) . The equation for the case of polarized electrons is: 



k^Vptl{ki,n,t) = h + h + h 



(67) 



where 



/i = ia 



i2ny 



E 



AYAZ 
R 



PZ'niZ) 



Pl,{k^)M".. 



a's[ (nYZ 



\nYZ 



::)-;n.,M;;f(:- 



(68) 



(27r) 



E 



AYAZ AY'AZ' Ak 



R 



R! 



6{k-k,)S{kcos0-kiCos0,)p-^,^{Z)p-;:,,AZ') 



2m:7^ ; "^^ 

' nYZ 



k " '^ \ n'Y'Z' 



:>M.,-M-(:- 



*■)«<■■;,(' "'^■'2' 

k '" " V n'Y'Z' 



.J^^^(^^ 






n'Y'Z' 
n'Y'Z' 



k) "'' \nYZ 



P%{h) 



(69) 



h = 



(27r 



)3^ 

n.7v 



AYAZ AY'AZ' Ak 



R 



R' 



S{R + k-R' - ki)6{Z + kcose- Z' -ki cos Oi) 



Sf p^JUz) - 6-;' pCiZ') 






nYZ 
n'Y'Z' 






n'Y'Z' 
nYZ 



fei 



Ps).\k^) 



(70) 
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'" " ^ nYZ 






fci 



P?, (fci 



P^' (fe) 



2m;,^ 



n'F'Z' 



fc 



, V",' n'Y'Z' 
, M „ ,} 
ki " ' \ nYZ 



fci 



pi! {k)p'ii{z) 5f -p%,„Az') 






h 






n'Y'Z' 
nYZ 



pl,„AZ') 5%, -pZ,^{Z) 






k / ''' \ n'Y'Z' 



fci 



p:, (fci) 



There are three terms in the rhs of equations (67). The first and the second terms describes redistribution only 
over polarization. The last one describes the general redistribution of the photons over quantum states (energies, 
momentum directions and polarization). 



2. Equation for the case of non-polarized electrons 

The kinetic equation in the case non-polarized electrons one can deduced from equation (67) by averaging over 
the spin states of the electrons. There is a relation between the distribution function of the electrons ,fn{Z) and the 
diagonal elements of the electron coherency matrix: 



UZ)=pIJZ)+pUZ), 
where p\^{Z) = p'^^{Z). The distribution function is normalized to the total number of the electrons: 



Y,J<iZUiZ) = N, 



Then equation for the case of non-polarized electrons is: 

kiYfill {ki ,ri,t) ^ J1+J2 + J3, 



where 



Ji 



{M- 



■E 



dYdZ fn(Z) 



pl^{k,)K€P 



J2 = 



i2ny 



E 



R 2 

dYdZ dY'dZ' dfc , 



nYZ 
nYZ 



p:nki)M::^ 



(nYZ 



R 



R' 



-5{k — ki)S{k cos 9 — ki cos 6*1) 



2Mo 



:''( 



nYZ 



\nYZ 



k ) "^1 



a 



E 



n'Y'Z' 

.ra"s'" (n'Y'Z' 

"^ \n'Y'Z' 
dYdZ dY'dZ' dfc 






\nYZ 

fjz)u,{z') 
2 2 

,,,", (n'Y'Z' 
k / '" " I n'Y'Z' 



(71) 
(72) 

(73) 
(74) 



h 



pl, (fci) 



kA^.s'fnYZ 
k) "'^ynYZ 



p'Mh) 



(75) 



•" (27r)3 '^JR R' k 

n.n' 

I [luz) fn'iz')] [r/;:f p:i(fei) + T;;!:;„p:;"(fci 

fniZ) 



5{R + k-R' - ki)5{Z + kcos0- Z' -ki cos 61) 



-fn'iZ') 



pl,{k) + 2TXyAk)Uz) 

T!„Ji(0.,{k^+T^lnpi[{k^ 



I'niZ') 



where we use a notation for the product of pairs of the scattering matrix elements: 



rpik —- A/f^^ 



n'Y'Z' 
nYZ 



k """[ n'Y'Z' 



(76) 



(77) 



After averaging over the spin states of the electrons the last term in the rhs of the kinetic equation simplifies. The 
(5-functions under the integrals in the rhs of equations (67) and (73) can be use to reduce the number of integrations 
and to simplify the sum. 
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3. Equation for the case of non-polarized rarefied electron gas 



Another form of the kinetic equations can be obtained in the case of rarefied electron gas. Neglecting in equation 
(73) the terms containing squares of the electron distribution function, we gets 



where 



Ki 



(27r) 



hVpll (fci, ri,t)^Ki+K2 + Ks, 



'^Y^l 



dYdZ fnjZ) 
R 2 






pl){k,)M:-l^ 



fnYZ 
\nYZ 



(78) 



(79) 



K2^0, 



(80) 



K, 



E 



dVdZ dY'dZ' dfc 



R 



R' 



S{R + k-R' - ki)d{Z + kcos0-Z' -ki cos6ii) 



'2t:c + Tf:(f6"{ki) + Tx,„p^"{k,) 



(27r)3 

n.n' 

xl{UZ)pi'{k) 

-frAz') [p%{h) (r,<;:f' +T,<;:fp:;'(fc)) +p^;"(fei) (t:'-,„ +T;i:;„p^;'(fc))] } . 

We notice that the first term of equation (73) has not changed, while the second term has now disappeared. 



(81) 



4- Kinetic equation in terms of two polarization modes 

In the case of non-polarized rarefied electron gas, equation (78) can be simplified further if one assumes the absence 
of correlations between the two linear polarization modes: p1(k) — p\{k) = 0. Then one can use only one polarization 
index for the diagonal elements of the coherency matrix: Pi{1i) = p\{k). The kinetic equation then get the form: 



k-^Vj)s^{ki,ri,t) = Li +L2 + L3, 



(82) 



where Li = and L2 = 0, and 

„2 



L, 



a 



{2^f 



E 



dydZ dY'dZ' dk 



R 



R' 



S{R + k-R' - ki)5{Z + kcose- Z' -ki cos 9i ) 



xTrsAfn(Z)pM[^+PsAkl)]-frAZ')psAkl)[l + Ps{k)]}. 



(83) 



This form of the equation is obvious and can be written immediately using physical arguments [22]. In this case 
the two modes are considered independently and a possibility of correlation between their phases is not taken into 
account. 



B. Kinetic equation in terms of Stokes parameters 



Equations (67) and (73) can be rewritten in terms of Stokes parameters. Transformation to this form can be 
done using trivial linear transformation. Elements of the coherency matrix {p){k)} and the Stokes vector TV = 
(ni, riQ, riu, fiv)"^ are connected by relations: 



ni^ipl+pl)/2, nQ^ip\-pl)/2, nu = {pj + pl)/2, ny = iipl - pl)/2, 
pl=ni + nQ, pl = ni-nQ, pl=nu-inv, pl=nu + inv. 



(84) 
(85) 
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1. General equation for the case of polarized electrons 



Using equation (67), one can find the kinetic equation in the ease of polarized electrons in terms of Stokes parameters: 



fci VTVi = /f + /f + /; 



3 ' 



where 



ir 



'(2 



'^rVJ 



dYdZ 
R 



f:,,XZ)TN,, 



(86) 



(87) 



/." = 



1 a' 



2 (27r)3 ^J R 



AYdZ AY'AZ' dfc 



R' 



6{k^ki)6{kcos9~kicosei)f^,^{Z)ff,',^,{Z') pN+J'"Ni , (88) 



p 1 a2 /• drdzdr'dz'dfc,, , , , X.. , . , , 

^3 = t;77;Z\^z2 I — 5 5^ -S{R + k- R' -ki)S{Z + k cos 9- Z' -ki cos Oi) 



2 (27r) 



n,n' 



R R' 



(89) 



X { [5frJL{Z) - S^J.'rZiZ')] UN, + 2/;;:(Z) [5-'" - /;,„„,(Z')] -R'N- 2f^,„„,{Z') [c" - /.<1(^)] n"N,}, 

where N = N{k) and Ni — N{ki) are the Stokes vectors, T, T' , T" , 7^, VJ and 7^" are 4x4 complex matrices 
acting like linear operators from the real 4-dimensional space to the real 4-dimensional space. The expressions for 
these matrices can be found in Appendix B. 



2. Equation for the case of non-polarized electrons 

The equation for the case of non-polarized electrons can be derived from equation (73): 

fciVATi = jP + jP + jP^ 
where 



'''-^ik^l 



AYAZUZ)^ 



(90) 



(91) 



Jr = ^ 



1 a^ 



2(27r) 



E 



dydZ AY'AZ' Ak 
~R 'R' k 



S{k — ki)d{kcos9 — ki cos 6*1) 



Uz)u{z') 



pN+P'Ni 



(92) 



^P 1 a2 ,-^ f AYAZ AY'AZ' Ak^. , „, , x..^ , ,, ^, , 
J3 = o77;:7;3 2^ / — T^ TT, -5{R + k - R' ~ ki)5{Z + kcos9 - Z' - fcicos^i) 



2(27r) 



R R' 



X \iuz) - f,Az')]nN, + uz) 



fn'iZ') 



TZ'N-fn'iZ'} 



UZ) 



n"Ni 



(93) 



where JF ^ JF' ^ JF" ^ TZ, TZ' and TZ" are 4x4 complex matrices that can be found in Appendix C. 

3. Equation for the case of non-polarized rarefied electron gas 

We derive the equation for the ease of rarefied electron gas by neglecting terms containing squared electron distri- 
bution functions in equation (90). The equation takes the following form: 



iiViVi = K[ + RP + RP, 



(94) 
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where 



^^-'J^I^J -H —^^^' (95) 



X|^ = 0, (96) 



p 1 a-' v- /"drdZdr'dZ'dfc,, , . . 



n.n' 



X [((/„(z) - fn'{z'))n - /„-(z')7^") ATi + /„(z)7^'A^] , (97) 

where J^, 7^, 7^' and 7?." are the same 4x4 complex matrices as in equations (90)-(93) and presented in Appendix C. 

V. SUMMARY 

We have deduced a kinetic equation for Compton scattering of polarized radiation in magnetic field. Polarizations 
of photons and spin states of electrons, the induced scattering and the Pauli exclusion principle were taken into 
account. The equations are written for both the coherency matrix and the Stokes parameters. Additional forms of 
the equations valid for the two polarization mode description of radiation is also derived. The equations for both 
polarized and non-polarized electrons were obtained. There are no significant (for the conditions in neutron stars 
atmospheres) limitations on the energies and the concentrations of the electrons and the photons. The assumptions 
made are usual for the kinetic theory and related to the typical time scales of the problem and do not limit significantly 
the applicability range. 

The equations describe the interaction of radiation and electrons in strong magnetic field up to about 10^^ G. There 
is no low limit on the B-field strength. At the same time, the derived equations become rather cumbersome in the 
case of weak magnetic field, because the electrons can occupy high Landau levels and there will be many terms in the 
rhs of the equation, where the summation over the Landau levels is carried out. On that other hand, in the case of 
strong magnetic field in neutron star atmospheres electrons typically occupy only ground Landau level (or only a few 
low levels), because of a rather low electron temperature and absence of high-energy photons. Therefore, the sums 
over n and n' has only a few terms, which simplifies the equations significantly. 

The most general form of the kinetic equation (67) has three terms. The second and the third terms there contain 
the products of two elements of the scattering matrix. These terms can be rewritten through the interaction cross 
sections. On the contrary, the first term in the rhs of equation (67) contains single elements of the scattering matrix. 
This term describes changes of the photon polarization with no corresponding changes in energy and the momentum 
direction. The polarization change term has the following form: the changes of the diagonal elements of the coherency 
matrix depend only on the non-diagonal elements and changes of the non-diagonal elements depend only on the 
diagonal ones. This term describes the rotation of the polarization plane when radiation can be well described only by 
the Stokes parameters, and it disappears, if the kinetic equation is reformulated in terms of two polarization modes. 
It is possible that this term provides correction to the depolarization in the region of vacuum resonance, which most 
likely will not be large because of a small optical depth of this region. 

The second term in the rhs of equation (67) describes redistribution of photons with only changes in polarization. 
It contains the products of the electron distribution functions and disappears in the equation (78) for rarefied electron 
gas. The last term in the rhs of equation (67) describes the general redistribution of photons over energy, directions and 
polarizations. This term can be simplified significantly for the cases of non-polarized electrons (73), rarefied electron 
gas (78) and for the two-polarization mode description of radiation (82). In the latter case this term is the only term 
in the rhs of the kinetic equation and coincides with previously known expressions. The derived equations form the 
basis for the construction of models of the radiative transfer in strongly magnetized neutron stars atmospheres and 
magnetospheres . 
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Appendix A: Transformation from density matrix kernel to distribution function for the electron in magnetic 

field 

The transformation from the kernel of density matrix to the distribution function for the case of charged particles 
in magnetic field is slightly more complicated than in the field-free case. We use the Landau gauge, with the electron 
momentum is described by two continuous components Y and Z (corresponding to the y- and z-coordinates) and 
one discrete component, which corresponds to the x-coordinate. A transformation from the coordinate representation 
to the momentum representation cannot be done using Fourier transforms. Instead one should use a transformation 
based on the Ritus eigenfunctions which diagonalize the mass operator of electrons in presence of the external magnetic 
field [32]. 

The kernel of the electron density matrix in the coordinate representation can be obtained from the kernel in the 
momentum representation: 



''"V. -) - s/ tStSi""'"'"""'"*"''"''*"""""' ' ("':- ' ■ '"'' 



where u = x — Y/h, u' = x' — Y' /b and <I'„(x) are the functions related to the parabolic cylinder functions: 

^n{x) = ^ D„{x). 

V27rv«' 
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These functions compose the complete system satisiying the relations 



/ $„(x)$„/(x)dx = (5;j', ^$„(x)$„(a;') =(5(a;' -, 



The inverse transformation from the kernel in coordinate representation to the kernel in momentum representation 
can be written in the following way: 



One can define the Wigncr function for the electrons in magnetic field: 

p%Jr, y,Z) = j di;^di;,e^(^'^«+^''^) / da;'$„(u')Pa' i^' , V - Vy/2, z - v,/2, x,y + Vy/2, z + «,/2), (A3) 

and rewrite it using the density matrix in momentum representation 






)y+(Z[-Zi)z]glY 



Yi+Y{ 

2 



.,(z_£i±£:)*.,,„„:y(;^'^:). ,a4) 

The inverse transformation from Wigner function to the kernel in momentum representation is 

p^'(r',r) = ^$„(y') /'dydZe-'[(^-«')^+(^-^')^lp-:(r,y,Z). (A5) 

From (A2) and (A5) one can obtain relation for kernel in momentum representation through the Wigner function: 

P-^'Cyz) = ^/^"'(g^{^^^| dVd3/e-(^'^'-^^+^'-'-^-)$„K^0^»(^)E^»"("0 (A6) 

To understand how to simplify equation (A6), let us consider a 1-dimensional problem, where the electron is 
described by its x-coordinate. The expansion of the density matrix kernel in coordinate representation is given by 
the following expression 

p{x',x) ^^<^n'{x')<i>n{x)pn,n'- (A7) 

The inverse transformation is 

Pn,n' = / dx'dx <^n'ix')<^nix)pix',x). (A8) 

Combining (A7) and (A8) one gets: 

pix',x) = E$„,(x')<I'„(x)ydx"'dx"$„.(a;"')<I'n(^")p(a;"',a;") ^J^I'nK^^OPnK^;), (A9) 

n.n' n' 

where 

p„,(a;) =^$„(x) /"dx"'dx"$„-(a;"')$„(x")p(x"',x") = j Ax"'dx"<^,,,{x"')6{x" - x)p{x"' ,x") = 

= f dx'^n'{x')p{x',x). (AlO) 
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If one considers the density matrix on a sufficiently short time scales, then there are no mixed states in the Landau 
levels and the coefficients Pn,n' in equation (A7) have only diagonal terms Pn,n' = ^n Pn,n- Therefore, 



p{x',x) = ^$„/(a;')$„(x)/9„,„. 



(All) 



Substituting this to equation (AlO), we get 

p„(x) = ^nix)pn,n- 

This proves that in a more general case p'^nr(''',Y', Z') = (f)„/(u)/9^„(Z'). This simplifies relation (A6): 



an 

ran 



d3,e^[(v'->').+(^'-^)^]$„(^)<j,„,(^)^^;^,(Z) 



YZ J (27r)2 

= ^/R;AZ^^/lUZ)SiY' - Y)5{Z' - Z)8ipi^,{Z') = RAZ)5{Y' - Y)5{Z' - Z)5i p^^Y, Z). 

This completes the proof of equation (26). 



(A12) 



(A13) 



Appendix B: Expressions for matrices in equations (87)— (89) 
The expressions for the matrix T — {Fij} in equation (87) can be written in the following form: 



-Fife — -Ffel — F; 



F23 = (^) - 
where we defined 



F.. 



AiD + iD), F, 
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m„i ^ 0, (fc, TO = 1,2,3,4), 

^23, ^34 -*((}>- (i», 



^4' 



-F.. 



O ^ K,' 



nYZ 
nYZ 



^^'3 ^ - 2ijlf 



F, 



43 



The elements of matrices T' ~ {-F/ } and T' ~ {Fl'A can be represented as 

Of"^,!! _L^,12 _|_^,21 I „,22n fPi2\ o/^,!! „,12 _l_ „,21 -,- ^,22\ 

^(7ii +721 ±7i2 ±722)> [f^J =2(7ii -721 ±712 T722)) 

^p,j _ ., ,11 , 721 ± 7?l ± 7I2), (^i:) = 2*(7n - 721 ± 7?l T 7I2), 

(t ) - 2 Q) [ill + ill ± 71^2^ ± 72^2^), (t ) = 2 (1) {^11 - ^11 ± ^11 T 72^^), 

(g) = 2 (1) (7?? + 7li^ ± 7^ ± 7l2^), (t ) = 2 (-1) (7?? - 7li^ ± 71^1 T 7l2^), 



-F: 



34, 



(Bl) 



(B2) 



(B3) 



where 





F" = 

^33 


„,ls /-Is „2s /• 

-7si - Gi - 7s2 - C 


2s 
s2 


(j- 1,2,3,4), 






Fi^^F^, = -jl!-C! + ili + Cl2, 


P" F" ^,ls /-2s „,2s aIs 
-^13 - -^31 - ^7s2 - Csl - Isl - Cs2 , 




Fi', = Fi[=^{Jli~a~J'/,+CU), 


F" F" ^,1« /'2s 1 ^,2s 1 /-Is 
-^23 - --^32 - -7s2 - Csl + 7sl + Cs2 , 




F!^, = -F^2 = ^illi Cst + Is! - C2), 


F34 = -Fi', = ^(-lU + a + Isi - Cs2) 


iL 


= ^-^ [nYZ 


fey """V"'^'^' 


k\ 

kij ' 


.31 -M^"'^Y"'^'^' 

•-fem a"k y^^/y'Z' 


kJ ""^ynYZ 



(B4) 



(B5) 



The matrix 7?. is a sum of the products of the elements of vector N (the Stokes parameters) and four matrices 

7?. = T^-ini + T^-qtiq + T^-unu + TZyny. (B6) 
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The elements of the matrices TZi^Q^uy are: 
(I 



— /ill _|_ a21 I a12 _|_ a22 



{R22^ 



- /l" in 421 I Al2 -^ a22 

— /i^-^ + /1]^2 + ^21 + ^221 



^31 

fl41 
fl42 



(I) 



(I) 



±cll + cll±cll + a 



12 

22' 



±Di\ + Cll±Dii + Cli., 



-R13 
-R14 

-R2; 
-R2. 



{!) 



(I) 






\R3i) 



Q)(±cll + c!^±c^ + ci 



Q) (TDll - ell T Dl 



c. 



22^ 
22 ^ 

22\ 
22 i: 



(B7) 



(Q) 
(Q) 
(Q) 
,(Q) 



_ /111 I /i21 _ /1I2 /.22 

R12 y — ^11 =■= ^12 ^21 + ^221 



^^21 V^^ _ 411 -p 421 _ ^12 , ^22 

i-R22i —^11+^12 ^21 ="=^22' 



ill 



22 



^^\\ + ^12 + "-^21 ^: 



'iV2 
22' 



(g^^) = ±i?21 ^ c-ll T ^2^2 - C-l^ 



-R13 

i?14 



(Q) 



(Q) 



a)(Sl\^TSl^^-i?2^?±B2^2^), 



(1) (±Cii + C?] T C]? - Q 



22 N 

22 ^ 



\R34J 

{rZ) ^""^ - a) (TDll ell ± D\l + Cll\ 



(B8) 



\Rv 



(U) 

- 4II + ^421 , ^12 , ^22 

_Ri2 y — ^21 =■= ^22 I ^11 =■= ^12 1 



\R22J 

(U) 



(R3i\ 

\R32J 



^21 + ^22 



.12 .22 

^11 + ^12i 



iU) 



-<12 



" ="=^21 + ^22 ^ "-^11 + ^12 

±Dll + Cll±Dll + CH., 



Rl3 
Rl4 

R2 
R2 

R33 
R34 



(U) 
(U) 
(U) 
(U) 



^Q){±All + All±All + All), 
-Q){BllTBll+BllTBll), 

{l){±cll + cil±cll + cll), 



{Z) -Q)iTDll-CllTDll-Cli), 



(B9) 



Rii 
Ri'. 

R21 
R22 

R31 
R32 



(V) 
(V) 
(V) 
,(V) 



— H ^21 T A22 + All ± A12) , 



H ^21 i ^22 + ^11 "F Ai 



22\ 

2J1 



H"F^2i ^ '-'22 i ^11 + ^ 



-(12n 



(£0 -*(Ti?2^^-C2\^±^??+Cl\^) 



(V-) 



f«23V^) 
1^24^ 



= Ci)(-^^^tA^}+A}^±A??), 

= Ci)(TS2\^+S2^2^±S??-Si^% 

= Ci)(-C,\iTCi+C|?±Cg), 



Z:)''^' = {l){Dll±Cll-DllTCli), 



(V) 



(BIO) 



where we introduced the fohowing combinations of the scattering matrices: 



Ai 



km 



<"/ 



n'Y'Z' 
nYZ 



M: 



a' I 



nYZ 
n'Y'Z' 



k\ ,,a, f n'Y'Z' 
uM^'AnYZ 



M: 



a" I 



nYZ 
n'Y'Z' 



(Bll) 



hm er" h 



n'Y'Z' 
nYZ 



ki\ ,, ( nYZ 
k) '^" I n'Y'Z' 



kj -'''[nYZ 



fc ^ ^^"''' 



nYZ 
"""^ ' n'Y'Z 



(B12) 



C: 



km 



K4 



n'Y'Z' 
nYZ 






M 



cr'/c 






'^'"(3-m) \ ^,Y'Z' 



(B13) 



£) 



km 
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nYZ 



k ™ I n'yz' 



(3-,) / n'F'^' 
^^^'^''^ ^ nYZ 



1 \ b^ct"; 



n'F'Z' 



(B14) 



The elements of the matrices it' are 



22 



I ij'^2 / ~ ^11 ^21 =■= ?12 =■= S22; 



«2l 






V«42y' 



<^21 

s;ii 



<r21 , ^12 <^ 
?12 ='=S21 +S; 



22' 



ell I <^22 , <^12 
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*11 _i_ ^22 -^ tl2 



e{2±^2iTe: 



22 1 



(£) 



(i)(«!5+«g±«;;±«g) 



«2 
«33 

«43 
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c22\ 
12 j 






(B15) 



where we defined 






m:,= 



n'Y'Z' 
nYZ 



MZ 



nYZ 
n'Y'Z' 



(B16) 



The elements of the matrices TZ" are 



where we defined 



"krr. 



m:j 





= 0l^ + el^±</)2^± 


& 


(«l) = a)(±0i^±pi^+<^?^+p^o 








[r'^J - [r-J ' {r'^J ~ 


-Q){Tc^l2±e'A + '^ll~ell), 






is A 

i'32j 




(R'is\ _ (-R'L\ 
' \R'L) - V «'i'J ' 
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'Z 


kA^^,(nYZ 
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fey V 
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(B17) 



(B18) 



Appendix C: Expressions for matrices in equations (91)— (93) 

Matrix J^ = {Fij} in equation (91) has the same form as in the case of polarized electrons (Bl), but where a' — a 
should be substituted in the definitions of (■[). Matrices P' and P" also have the same forms as in the case of polarized 
electrons given by equations (B3) and (B4), where a' ^ a and a'" = a" should be used in the expression for 7*Jj^ and 

/■ik 
^jm ■ 

Matrix it can be represented by expression (B6). The elements of TZj can be represented through tensors T*,^ 
defined by relation (77): 
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The elements of TZq can be represented as 
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The elements of TZu are 



u 

{rA ) '""^ = ^n ± Ti'l + T^i ± T|2^ , («;: ) ^""^ - (!) (^iT ± T^2 + Til ± ni 



rW-R'^\ Rr.'-o, 4T-R'l(\ RiV=o, 

R^lP=^{T^^+nl)^ <)=z(ri| + r,i,i), 4^)=0, <) = T^f + r,^| + T,Y + Tf^^ . (C3) 

The elements of TZv are 

-RllV _ T^12 I rp22 rpll _^ rp21 /-RlsA "^ _ /'lW_l_rp22 , rpl2_rp21 rpUN 

_Ri2j — -'11 =■= Ji2 — -t2i + ^22 1 \^fli4j — \i) V^-'ll + -'12 + -^21 ^ -'22 j 7 



-R22y ~ V«liy ' \R'2i) ~ ^' 



•-31 " -"-13 ' ^32 ~ '-'' ^33 ~ ^11 ' ^34 ~ 

y{V) _ rp22 , rTn21 n(V) _ rTnl2 , rpH tj{V) _ „ jy{U) _ ■ ( rpl2 . rp22 rpll rp2l\ 



-n-41 — ^-'ll+-'21: -^^42 —^-'12 +-'22: -'^43 ^ ^' -^44 — ^ 1,-' 11 + -^ 12 ^ -^ 21 ~ -* 22 j • \'^^) 

The expressions for the elements of TZ' are 

(^1:) = Til + Tli ± Tif ± Til (§:) = Q) {Tli + Tii ± Til ± TH) , 

(R'2l\ _ rpll _ rp21 1 rTnl2 _ rTn22 ^■^23 A _ rpl2 _ rp22 _1_ rpll rTi21 

V^22/ ~ " -^12 =t ^21 + -^22 1 Vfl24/ ~ " -^12 =t ^21 + -'22 7 

R'3l\ _ rp2l I rpll I rp22 _1_ rpl2 /^-RsS A _ f l\ /'Tn22 , rpl2 , rp21 , rpll 

-^11 + -^12 =■= ^21 =■= J22 ; Ifl' i — uM-'ll + -^12 =■= i21 ="=1; 



\^32/ 



22, 



(S:) - ^ (^n - ^1^2^ ± Til T Tii) , (^ = (1) {±Tii T TH + TH TH) , (C5) 

and the elements of TZ" are as follows: 

K>" — K>" R" — R" — R" — R" — n R" — R" R" — R" R" — R" ('r'R^ 



